Flow boundary conditions for chain-end adsorbing polymer blends 
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Using the phenol-terminated polycarbonate blend as an example, we demonstrate that the hydro- 
dynamic boundary conditions for a flow of an adsorbing polymer melt are extremely sensitive to the 
structure of the epitaxial layer. Under shear, the adsorbed parts (chain ends) of the polymer melt 
move along the equipotential lines of the surface potential whereas the adsorbed additives serve as 
the surface defects. In response to the increase of the number of the adsorbed additives the surface 
layer becomes thinner and solidifies. This results in a gradual transition from the slip to the no-slip 
boundary condition for the melt flow, with a non-monotonic dependence of the slip length on the 
surface concentration of the adsorbed ends. 

PACS numbers: 83.80.Sq, 61.20.Ja, 47.27.Lx 
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I. INTRODUCTION 

The equations of continuum fluid mechanics are in- 
complete without appropriate boundary conditions. In 
most situations it is required that both normal and tan- 
gential components of the relative fluid velocity vanish 
at the surface, 1 This, so called stick or no-slip bound- 
ary condition, has successfully accounted for most of the 
experimental facts. It is, however, empirical by nature: 
there are no theoretical arguments in favor of the no- 
slip; moreover, it has been known since Maxwell's times 2 
that even a simple kinetic theory of gasses predicts the 
non-zero value of the tangential velocity at the wall. 

Providing some of the insight into the question, kinetic 
theory fails already for simple fluids adjacent to a rigid 
solid. In a more general context, multiple scattering from 
the individual wall molecules remains the major problem 
of most analytical theories. In this situation computer 
simulation techniques are able to advance our knowledge 
of the processes occurring at the surfaces. 

Indeed, molecular dynamics (MD) simulations demon- 
strated that, in case of simple fluids, the flow boundary 
conditions are sensitive to the fluid epitaxial ordes^ as 
well as the wall structure** The situation, however, be- 
comes much more complicated for polymers at the sur- 
faces, because of the much richer molecular arrangement 
(e. g. formation of brushes^ various adsorbed and de- 
pleted layers^) as well as much stronger correlation be- 
tween the atoms right at the wall with the rest of the 
surface layer i£*L 

Recent MD studies of the end-adsorbing polymer 
melts, performed with a novel quantum based multiscale 
approach for the surface/polymer interaction, showed 
that at least two mechanisms contribute to the hydro- 
dynamic boundary conditions— The attached parts of 
the chains scatter on the surface potential while moving 
along its equipotential lines. This induces the density 
and the chain conformation modulation in the adsorbed 
layer, and energy is lost from these modulations through 
the coupling to the thermostat, similar to the situation 
observed for a generic model of adsorbed surface layersii*> 
On the other hand, single-end grafted chains of polymer 



brushes undergo a coil/stretch transition and disentan- 
gle from the melt at a given shear rateji2ii2ii»ii£i favoring 
the slip boundary condition. Though investigated for the 
special case of end adsorbing polycarbonate on nickel, a 
similar scenario will occur for block copolymers with an 
adsorbing block. 



Valid for monodispersed melts, this picture does not ac- 
count for the usual melt polydispersity, which is an out- 
come of all synthetic polymerization reactions. Moreover, 
in many situations the self-blending of polymer melts is 
used— A small amount of a lower weight polymer im- 
proves the melt processability without significantly af- 
fecting its mechanical properties. It, however, alters the 
structure of the surface layer— modifying the hydrody- 
namic boundary conditions for the melt flow. These 
changes shall be taken into account at the later stages 
of melt processing 



In this work we focus on the hydrodynamic bound- 
ary conditions for polymer blends adsorbed on a solid 
substrate. For this purpose we consider a particular sys- 
tem, that is bisphenol-A polycarbonate (BPA-PC) melt 
sheared over a (111) nickel surface. Our choice of poly- 
carbonate as a test system is twofold. First, it is widely 
used in various applications^ and, therefore, has been 
intensively studied both experimentally 20 and theoreti- 
ca lly ( 2i 4 22 i 23 J 24 i 25 j n particular, shear of monodispersed 
melts' and the structure of the static epitaxial layer for 
monodispersed melts^SiSi and blends*£ have already been 
considered. Second, although it is a specific system, it 
has a number of important generic features, which also 
apply to a realistic description of many polymer brushes 
composed of block copolymers. 



Our prime goal is to relate the structure of the ad- 
sorbed layer, which changes in the presence of the low 
molecular weight additive, to the hydrodynamic bound- 
ary conditions for the melt flow, specified by the slip 
length and the surface friction coefficient. 
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FIG. 1: Chemical structures: (a) phenol-terminated 
bisphenohA polycarbonate (in this study n = 1,5,20); (b) 
diphenyl carbonate (DPC); (c) phenol. 



II. SIMULATION DETAILS 

We consider four types of polymer mixtures. The host 
polymer (major component) is the phenol terminated 
BPA-PC of N% = 20 repeat units. The second (minor) 
component is one of the following: BPA-PC of N2 = 5 
repeat units, one repeat unit, diphenyl carbonate (DPC), 
or phenol, all shown in Fig. ^ For all mixtures, we use 
n\ — 400 chains of the major component; the number of 
molecules of the minor component, ri2, is then adjusted 
to provide (approximately) 5% of the total weight of the 
system. Exact numbers are given in Table [I] 





ni, ni 


L z ,a 


20 


400 


80.91 


20,5 


400, 76 


85.15 


20,1 


400, 239 


84.97 


20, DPC 


400, 521 


84.69 


20, Phenol 


400, 1186 


83.79 



TABLE I: Studied systems, a = 4.41A 

To simulate the melt, we use the previously devel- 
oped coarse-graining model, in which each monomer is 
replaced by four beads that correspond to isopropylidene, 
carbonate, and the two linking phenylenes. Interaction 
potentials, bead sizes and coarse- graining procedure are 
described in Refsi^&^l. The bead-wall interaction poten- 
tials are obtained from ab initio density functional calcu- 
lations. All internal beads experience strong increasing 
repulsion below 3.2 A, either due to the nature of the 
molecular interaction or due to the steric hindrance by 
the other beads. Only the chain ends absorb on the wall. 
The wall-end interaction potential is expanded in 2D re- 
ciprocal lattice space of (111) nickel surface and has the 
following form 



U(x,y,z) = ^2lli(z) fi(x,y), 
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The interaction potential obtained using ab initio cal- 
culations is well reproduced by the following set of pa- 
rameters: e r = 1.5eV, e Q — 0.7eV, e\ = — 7/45eV, and 
e 2 = -2/45 eV. For details, see Refs. 10 i 24 i 28 . 

The melts are confined to a slit pore of thickness L z 
with the walls perpendicular to the z axis. Periodic 
boundary conditions are employed in x and y directions. 
The x and y box dimensions are set to L x = 22.23 a and 
L y = 21.72cr with a = 4.41 A, which corresponds to a 
(111) hexagonal lattice of nickel with 39 and 22 unit cells. 
The number density of beads is n = 0.85, which corre- 
sponds to 1.05 g/cc, the experimental density at the pro- 
cessing temperature, 570 K. The units are chosen such 
that k B T = 1 with T = 570 K and a is unity^ 2 ! 

Starting configurations are generated by randomly 
placing the chains in the simulation box. A short run 
is then used to remove the bead-bead overlaps^ The 
production run is performed in the NVT ensemble with 
Langevin thermostat with friction 0.5t -1 , where r is 
the unit of time in the simulations. The thermostat is 
switched off in the shear direction. The velocity- Verlet 
algorithm with the timestep At — 0.005r is used to in- 
tegrate the equations of motion. After equilibration for 
about 2 x 10 5 t, the shear is applied by moving the top 
and bottom walls in opposite directions at a constant ve- 
locity v Wl so that the shear rate is 7 = 2v w /L z . The wall 
velocity is the same as in our previous studies of monodis- 
persed BPA-PC meltsjiS v w r/a = 0.01. The correspond- 
ing shear rate can be obtained from the time mapping 
It 25 ps, see Ref. 2 — for details. Using this mapping we 
obtain 7 10 7 s _1 . Taking into account that the average 
chain length in a BPA-PC melt is N ~ 70 the corre- 
sponding chain reptation time Td ~ N 3A is almost by 
two orders of magnitude larger than that of the N = 20 
chains. Equivalently, the shear rate would be reduced to 
10 5 s^ 1 , close to the value used for industrial processing 



of BPA-PC, 



10 4 



Note that smaller shear rates 



(1) 



either have no significant effect on the adsorbed layer or 
are difficult to analyze, due to significant error bars for 
the velocity profiles. 



III. RESULTS 



where i = 0,1,2 corresponds to the reciprocal vectors of 
different lengths, f = 1, /1 = cos(x — y) + cos(x + y) + 
cos 2y, and fa = cos(a; — 3y) + cos(x + 3y) + cos 2x, where 



A. Blends with different additives 

We first have a look at the structure of the surface 
layer. Fig.|2shows a typical chain end density profile for a 
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monodispersed melt (N = 20) . The chain end density has 
a sharp peak next to the wall, which is due to the strong 
adsorption of the ends; the region with no ends follows; 
finally, the bulk concentration is reached at a distance of 
the order of the radius of gyration, R g w Q.2a = 27 A. 
The inset illustrates how this profile changes in the pres- 
ence of the additives: the shorter the chains of the minor 
component, the more of them adsorb on the wall. This 
of course results in a decrease of the number of the ad- 
sorbed chain ends of the major component, in agreement 
with our earlier studies^ except for the phenol additive, 
for which we previously observed a weak increase in the 
number of the adsorbed chain ends of the major compo- 
nent compared to the DPC casei^ 





20:1 

Mixture 



20:DPC 20:Phenol 



z/a 



25 



FIG. 2: Typical chain end density profile (N = 20). The 
inset depicts the change in the surface density (number of the 
adsorbed chain ends per unit area) for the minor and major 
components. The log scale is used to show the bulk density 
and the depleted region on the same plot. 

The representative snapshots of the systems with and 
without shear are shown in Fig. [3J Snapshots (a) and 
(c) illustrate the thinning of the adsorbed layer due to 
the decrease in the number of the adsorbed chain ends of 
the major component^ After shear is applied, the one- 
end attached chains disentangle from the melt, stretch 
and form a thin lubricating layer between the bulk and 
the chains adsorbed with two ends (see snapshots (b) 
and (d)). Corresponding changes can also be detected 
by observing the center of mass profiles^ 

Having resolved the structure of the adsorbed layer, we 
shall turn our attention to the velocity profiles, which are 
shown in Fig. 0] For all mixtures, the profiles have fea- 
tures similar to those of monodispersed systems^ Next 
to the wall the velocity is practically constant, since the 
adsorbed chains are dragged by the wall. Immediately 
after the plateau, the velocity profile becomes a linear 
function of z. The velocity of the beads at the wall, 
v s , is smaller than the wall velocity, v w , i. e. the chain 
ends slide over the wall, moving between the hollow and 



FIG. 3: Snapshots of the mixtures: 20:5 without (a) and 
with (b) shear, 20:phenol without (c) and with (d) shear. 
Polymer chains are divided into three populations: chains 
which adsorb both ends (green), only one end (red), and no 
ends on the surface (blue). 



bridge sites of the potential (the difference in the adsorp- 
tion energies of these two sites is rather small, of the 
order of 2kT, see the contour plot in Fig.[5J). In case of 
the monodispersed system and 20 : 5 mixture we have 
v s /v w ~ 0.5. This ratio increases for the 20 : 1 mixture 
and reaches 1 for the diphenyl carbonate and phenol ad- 
di" 




FIG. 4: Velocity profiles normalized to the wall velocity. 
Only a part of the profiles in the vicinity of the wall is shown. 
The inverse of the surface friction coefficient /3 _1 ~ l—v s /v w 
(circles) and the shear stress f x /S (squares) are shown in the 
inset. To calculate the velocity profiles we average over all 
beads, independent of the nature of a particular bead. 

We can write the stationary shear stress t xz (force per 
unit area) as^ii 

T X z = fx/S = j3{v w -v s ), (2) 

where (3 is the friction coefficient between the adsorbed 
layer and the wall, S is the surface area. Since the den- 
sity of the minor component in the bulk is small (5% 
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in weight and most of it is adsorbed on the walls) the 
bulk properties of the melt do not change significantly 
from mixture to mixture. Indeed, the shear stress t xz 
is practically independent of the molecular weight of the 
additive (see the inset to Fig.^J. Therefore, the change of 
v s is due to the different values of the friction coefficient 
/3, which increases significantly for short additives. 

This increase is, in fact, rather unexpected: the ad- 
sorbed additives reduce the number of the attractive 
sites available for the major component and, in principle, 
shall screen the effective interaction of the melt with the 
surface. Weaker interaction should result in a smaller, 
compared to the monodispersed melt, friction coefficient. 
This, however, does not happen in practice: the systems 
with the phenol or diphenyl carbonate additives, which 
basically cover about 80% of the wall, have the biggest 
friction coefficient. 




FIG. 5: Snapshots of the chain ends adsorbed on the walls. 
Blue beads: major component, yellow beads: minor compo- 
nent. The bead sizes correspond to their van der Waals radii. 
The triangular grid depicts the lattice of the (111) nickel sur- 
face. Contour plot illustrates the bead-wall potential used in 
simulation, together with the possible paths for the adsorbed 
beads. The contour levels are in electron volts, 1 eV ~ 20 kT. 

To understand why this happens, let us have a look at 
the structure of the beads adsorbed on the surface. Fig.OJ 
shows the snapshots of the representative mixtures. The 
difference between the systems is now clear: the adsorbed 
layer of the 20:5 mixture has a structure similar to a two- 
dimensional (2D) gas; the packing of the 20:1 mixture is 
more dense, but still disordered, similar to a 2D liquid; 
for the 20:DPC as well as 20:Phenol (not shown) mixtures 
the concentration of the adsorbed chain ends is so high 
that they form a 2D crystalline layer. The hexagonal 
lattice of this layer has the lattice constant of 2a, where 
a is a lattice constant of the (111) nickel surface, which is 



of the order of the van de Waals diameter of the adsorbed 
beads, i. e. the close-packed surface layer is not frustrated 
energetically, or, in other words, the surface potential and 
the surface layer have commensurable lattices. 



B. Different concentrations of the phenol additive 

To quantify our results even further, we have also stud- 
ied mixtures with different concentrations of the phenol 
additive. In these mixtures the surface densities of the 
adsorbed chain ends and phenol molecules are monotonic 
functions of the total number of phenol molecules in the 
system. Hence, they can be varied with a good precision, 
which helps to analyze the systems in a more systematic 
way. 

The systems were prepared in a similar manner: the 
number of the phenol molecules was adjusted to provide 
different percentages of the total weight of the system, 
in the range from 1% to 20%. Because of the finite size 
of the simulation box and phenol adsorption, the con- 
centration of the phenol molecules in the bulk is always 
smaller than the one used for the system preparation. 
We therefore used this number (from 1% to 20%) only to 
label a particular system. The important quantities are 
of course the surface density of the adsorbed chain ends 
and phenol molecules. 

In addition to the velocity profiles we have calculated 
the two-dimensional radial distribution function of the 
adsorbed molecules 

v ' \ i=Q j=i+l I 



where r i3 = y/(xi - xj) 2 + (yi - y-j) 2 , S{r) = n[(r + 
Ar) 2 — Trr 2 ] is the bin area, Ar = r max / N^ins is the bin 
width, pij (r) = 1 if r < < r + Ar and zero otherwise. 
The sum is performed over the adsorbed beads (chain 
ends and phenol molecules) only, i. e. Zij — z wa ii < cr; n 
is the number of the adsorbed beads; < ... > denotes the 
ensemble average. 

To further characterize the quality of the hexagonal 
packing we have also calculated the orientational order 
parameter 31 

Se = ( ' ^ n ) (4) 

where (f>ij is the angle between the vector r^ and a given 
axis in the xy plane, </y = 1 if r^ belongs to the first peak 
of the radial distribution function (in our case between 
0.5ct and 1.5<r), and zero otherwise. Note that S& = 1 
in the case of a perfect hexagonal order whereas = 
indicates the complete lack of such order. 

The radial distribution functions are shown in Fig. |SJ 
The increase of the total number of the adsorbed 
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FIG. 6: Two-dimensional radial distribution function of the 
adsorbed chain ends. Inset shows the surface densities (num- 
ber of adsorbed beads per unit area) of the adsorbed ends as 
a function of the phenol concentration. 



molecules results in a gradual solidification of the sur- 
face layer. The 4% system (p s — 0.775er~ 2 ) has a two- 
dimensional solid layer formed at the surface, as it can 
be seen from the peaks at 21 and where I is the posi- 
tion of the first peak. Further increase of the amount of 
adsorbed molecules merely improves the hexagonal close- 
packing and, after the 5% system (p s = 0.851a -2 ) the 
close packing of the surface is reached, i. e. the radial 
distribution function does not change anymore. The ra- 
tio between the number of the adsorbed chain ends and 
the adsorbed phenol molecules, however, still decreases, 
as it can be seen from the inset to Fig. [§] In case of the 
blends, the solidification of the surface layer occurs for 
the 20:DPC mixture. 

The force on the wall and the velocity of the adsorbed 
layer are shown in Fig. Q together with the orientational 
order parameter Sq. The orientational order parameter 
Sg confirms that the solidification of the surface layer oc- 
curs for the 4% system, and the saturation for the 5% 
system, in agreement with the behavior of the radial dis- 
tribution functions. 

From the dependence of the relative velocity of the 
adsorbed layer on the surface density of the adsorbed 
molecules one can see that once the epitaxial layer is in 
a solid state, it follows the wall, i. e. we have the stick 
boundary conditions for the surface layer. However, be- 
fore the solidification, the velocity of the adsorbed layer 
is smaller than the wall velocity and is roughly propor- 
tional to the total concentration of the adsorbed chain 
ends. 

Similar to the relative velocity, the friction force on the 
wall increases with the increase of the total amount of the 
adsorbed beads. However, it does not reach a constant 
value at the point of solidification of the surface layer, as 
the relative velocity does, but starts to decrease, with the 
decrease rate proportional to the ratio of the number of 
the adsorbed chain ends to the number of the adsorbed 



phenol molecules. This is due to the fact that the amount 
of the long chains adsorbed on the wall still decreases and 
the surface layer disentangles from the bulk of the melt. 

Finally, we have also calculated the slip length, which is 
defined as the extrapolation length of the linear velocity 
profile in the bulk to zero and is often taken as a param- 
eter in macroscopic descriptions of a fluid flow. The con- 
centration dependence of the slip length is shown in the 
inset to Fig. [7| Two regimes can be clearly seen. First, 
at a low concentration of the adsorbed chain ends (e. g. 
in case of a pure melt, N = 20) the melt slips over the 
wall. If we increase the surface concentration of the ad- 
sorbed beads, the slip length decreases and becomes zero 
at some concentration, i. e. we effectively have the no-slip 
boundary conditions for the melt flow. This point, how- 
ever, does not coincide with the actual locking of the sur- 
face layer, which we observe only when it solidifies. The 
reason can be seen from the shape of the velocity profiles: 
let us denote the thickness of the adsorbed layer as 6 and 
the velocity of this layer as v s . Then the bulk velocity can 
be written as v = v s (L z — 2z)/(L z — 25). The slip length 
b, obtained from the condition v(z = —L z /2 — b) = v w , 
reads 

b=(v w /v s -l)L z /2-Sv w /v s . (5) 

It can be seen that two mechanisms contribute to the to- 
tal slip length b. The first one is due to the apparent slip 
of the adsorbed layer over the surface. The second, neg- 
ative, contribution is due to the finite thickness of the 
adsorbed polymer layer. As we add the additives, the 
adsorbed layer shrinks, giving rise to a higher slip; the 
relative velocity of the adsorbed layer, however, drops 
down much faster. Combined, these two mechanisms re- 
sult in a zero slip for some intermediate concentration of 
the adsorbed chain ends and a negative slip for higher 
surface concentrations. 

After the solidification, the adsorbed layer follows the 
wall, i. e. v w = v s . However, if we further increase 
the amount of the phenol additive, the ratio between the 
number of the adsorbed chain ends of the long chains and 
phenol molecules will decrease, as well as the thickness of 
the adsorbed layer. As a result, the (negative) slip length 
will start to increase again. 

IV. DISCUSSION AND CONCLUSIONS 

Let us first turn to a brief discussion on the possible 
mechanisms of friction in our system. Recalling that the 
adsorbed chain ends move along the equipotential lines 
of the surface potential, we conclude that the role of the 
adsorbed additives is twofold: apart from screening the 
interaction of the melt with the wall, they also serve as 
additional obstacles which block the possible paths for 
the chain ends of the major component. If wc ignore 
weak entanglement of the chains in the adsorbed layer as 
well as the motion of the adsorbed chain ends of the mi- 
nor component over the surface, the situation reduces to 
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FIG. 7: Normalized (to the wall velocity) velocity of the ad- 
sorbed layer v s /v w (circles), orientational order parameter Se 
(squares), and the force per unit area on the wall (triangles). 
Inset: slip length, defined as a distance at which the melt 
velocity extrapolates to zero. 

a two-dimensional site percolation problem. Each vertex 
is designated open or closed at random, with probability 
p to be closed which is proportional to the concentration 
of the chain ends of the minor component. Percolation 
theory predicts that for the site percolation on a hexago- 
nal lattice the critical point of the percolation probability 
is p c = 1/2. 32 Assuming that the transition happens for 
the 20 : 1 system (see the inset of Fig. we obtain 
p c » 0.47. Of course, due to the applied shear we have 
an oriented or directed percolation: our lattice sites shall 
be assigned particular orientations, along which the per- 
colative paths are biased. 

On the other hand, gradual solidification of the ad- 
sorbed layer results in the increase of the friction co- 
efficient, due to the collective motion of the adsorbed 
molecules: once the surface layer is solidified, the ad- 
sorbed molecules cannot move in the direction perpen- 
dicular to the shear direction. However, hopping between 
the hole and bridge sites always involves a motion per- 
pendicular to the shear direction, or, in other words, the 
hexagonal symmetry of the lattice forbids the motion of 
the beads along the shear. 

Once the surface layer is in a solid state, the relative 
velocity of the adsorbed layer does not change anymore. 
However, the ratio between the adsorbed chain ends of 



long and short molecules still changes. This affects the 
thickness of the adsorbed layer as well as its entangle- 
ment with the rest of the melt. As a result, the shear 
force on the wall decreases with the increase of the con- 
centration of the additive. Similar scenario has already 
been discussed in the framework of the mean-field the- 

Qr y 12.13.14.15 

To summarize, we studied the effective boundary con- 
ditions for a polymer blend adsorbed on a structured 
surface. The slip boundary condition observed for a 
monodispersed melt changes to the no-slip at some con- 
centration of the additive. Further increase of the con- 
centration of the additive at the surface leads to the so- 
lidification and locking of the motion of the adsorbed 
surface layer. 

Finally, we would like to comment on the importance 
of the multiscale modeling methods we employ in our 
studies. As underlined in the introduction, we consider 
a particular system, which is of high relevance to dif- 
ferent fields of modern technology. Our results, despite 
the particular system considered, are the direct expres- 
sion of the interplay between specific (electronic based) 
molecule-surface interactions and global statistical and 
dynamical properties of the system, i. e. the interplay 
between the different scales is the crucial ingredient of 
the description. The electronic and molecular resolution, 
implicit into the parameterization of the model, allows 
for a level of analysis which is beyond any other exist- 
ing coarse-grained models or mean-field approaches. In 
this sense, this study goes beyond the specificity of the 
system considered and calls for extensions to other sys- 
tems, and experimental tests. The present approach is 
a route which allows to detect important, otherwise not 
accessible, properties and sets a new link between theory, 
experiments and technology. 
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